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Statistical theory of shot noise in quasi-lD Field Effect Transistors in the presence of 

electron-electron interaction 
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We present an expression for the shot noise power spectral density in quasi-one dimensional 
conductors electrostatically controlled by a gate electrode, that includes the effects of Coulomb in- 
teraction and of Pauli exclusion among charge carriers. In this sense, our expression extends the 
well known Landauer-Biittiker noise formula to include the effect of Coulomb interaction inducing 
fluctuations of the potential in the device region. Our approach is based on evaluating the sta- 
tistical properties of the scattering matrix and on a second-quantization many-body description. 
From a quantitative point of view, statistical properties are obtained by means of Monte Carlo 
simulations on an ensemble of different configurations of injected states, requiring the solution of 
the Poisson-Schrodinger equation on a three-dimensional grid, with the non-equilibrium Green's 
functions formalism. In a series of examples, we show that failure to consider the effects of Coulomb 
interaction on noise leads to a gross overestimation of the noise spectrum of quasi-one dimensional 
devices. 

PACS numbers: 73.50.Td, 73.63.Nm 



I. INTRODUCTION 

As quasi one-dimensional field-effect transistors 
(FETs), based for example on Carbon Nanotubes 
(CNTs) or Silicon Nano Wires (SNWs), are increas- 
ingly investigated as a possible replacement for conven- 
tional planar FETs, it is important to achieve complete 
understanding of the properties of shot noise of one- 
dimensional conductors electrostatically controlled by a 
third (gate) electrode. Shot noise is particularly sensitive 
to carrier-carrier interaction, that in turn can be particu- 
larly significant in one-dimensional nanoscale conductors, 
where electrons are few and screening is limited 1]. 

Low frequency 1// noise in quasi one-dimensional con- 
ductors has been the subject of interest for several au- 
thors 043; whereas few experimental papers on shot 
noise have recently been published fs', ^ . 

Due to the small amount of mobile charge in nanoscale 
one- dimensional FETs, even in strong inversion, drain 
current fluctuations can heavily affect device electrical 
behavior. Of course, noise is an unavoidable and unde- 
sirable feature of electron devices, and its effect must be 
minimized or kept within tolerable levels for the oper- 
ation of electronic circuits. From a more fundamental 
point of view, it is also a rich source of information on 
electron-electron interaction, which cannot be obtained 
from DC or AC electrical characteristics. 

The main sources of noise are injection from the con- 
tacts into the device region, through the random occu- 
pation of states around the Fermi energy at the contacts, 
and partial transmission of electrons through the con- 
ductor, which gives rise to the so called partition noise. 
The main types of interaction that have a clear effect on 
noise are Pauli exclusion, which reduces fluctuations of 
the rate of injected electrons by limiting the occupancy 
of injected states, and Coulomb repulsions among elec- 
trons, which is the cause of fluctuations of the potential 



in the device region, that often suppress, but sometimes 
enhance the effect of fluctuations in the rate of injected 
electrons. 

The combined effect of Pauli exclusion and Coulomb 
repulsion on shot noise has been investigated in the case 
of baUistic double gate MOSFETs 0, in planar MOS- 
FETs 8] and in resonant tunneling diodes |9l-[ll|. There 
are still few attempts to a complete quantitative un- 
derstanding of shot noise in ballistic CNT- and SNW- 
FETs. Indeed, when addressing a resonant tunneling 
diode one can usually adopt an approach that exploits 
the fact that the two opaque barriers break the device in 
three loosely coupled regions (the two contacts and the 
well), among which transitions can described by Fermi 
golden rule, as has been done in Refs. QrJJJ. This is 
not possible in the case of a transistor, where coupling 
between the channel and the contacts is very good. 

Another important issue is represented by the fact 
that the widely known Landauer-Biittiker's noise for- 
mula [ij, does not take into account the effect of 
Coulomb interaction on shot noise through potential fluc- 
tuations. Indeed, recent experiments on shot noise in 
CNT-based Fabry Perot interferometers show that 
in some bias conditions many-body corrections might 
be needed to explain the observed noise suppression. 
Other experiments show that at low temperature sus- 
pended ropes of single-wall carbon nanotubes of length 
0.4 pLTn exhibit a significant suppression of current fluc- 
tuations by a factor smaller than 1/100 compared to full 
shot noise However, this experimental result is not 
supported by a convincing interpretation, since possible 
explanations extend from ballistic transport in a small 
number of tubes within a rope, to diffusive transport in 
a substantial fraction of the CNTs. 

In this work, we present an expression for the shot 
noise power spectral density of ballistic quasi-one dimen- 
sional channels based on a statistical approach relying on 
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quantities obtained from Monte Carlo (MC) simulations 
of randomly injected electrons from the reservoirs. The 
expression is derived within the second quantization for- 
malism, and simulations are based on the self-consistent 
solution of the 3D Poisson and Schrodinger equations, 
within the non-equilibrium Green's function (NEGF) for- 
malism jisl . 

Our proposed expression generalizes the Landauer- 
Biittiker's noise formula including the effects of Coulomb 
interaction, that is significant for a large class of devices, 
and in particular for one-dimensional conductors. 
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II. THEORY 

According to Milatz's theorem [IB], the power spec- 
tral density of the noise current in the zero frequency 
limit can be written as S{0) = \mii/^_^^ [2/i^ • var{I)], 
where v is the injection rate of a carrier from a contact 
and var{I) is the variance of the current. According to 
Ref . [l7| ■ v can be expressed as = AE/{2nh) where AE 
is the energy discretization step, i.e. the minimum energy 
separation between injected states. Indeed, the contribu- 
tion to the current of a transverse mode in the energy in- 
terval AE can be expressed in the zero temperature limit 
by the Landauer-Biittiker formula as (dl) = e/ (27r/i) AE. 
On the other hand {dl) = ev, from which v = AE/ {2nh) 
derives. Finally, the power spectral density of shot noise 
at zero frequency can be expressed as: 

SiO) ^ lim - var(I) = lim AnfC^^^^ . (1) 

The variance of the current can be derived by means of 
the second quantization formalism, which allows a con- 
cise treatment of the many-electron problem. 

Let us consider a mesoscopic conductor connected to 
two reservoirs [source (S) and drain (D)], where electron 
states are populated according to their Fermi occupa- 
tion factors (Fig. [T]). For simplicity, we assume that the 
conductor is sufficiently short as to completely neglect 
inelastic scattering events. Thermalization occurs only 
in the reservoirs. At zero magnetic field and far from the 
interacting channel, the time-dependent current operator 
at the source can be expressed as the difference between 
the occupation number of carriers moving inward (Ng^) 
and outward {Ng^^) in each quantum channel m jisj : 

= ^ E / ^ [NLiE, t) - NgjE, t)] , (2) 

where 

^5™ (E, t)^ j d{nw)b+^ {E) bsm {E + t^) e--* . (3) 

The introduced operators a^smi-E) and asm{E) create 
and annihilate, respectively, incident electrons in the 



FIG. 1; Annihilation operators for ingoing {as ) and 

outgoing electron states {bsm, bon) in a two terminal scatter- 
ing problem (m — 1, Ws; n — 1, Wd)- 

source lead with total energy E in the channel m (Fig.[T]). 
In the same way, the creation bg,j^{E) and annihilation 
i>Sm{E) operators refer to electrons in the source contact 
for outgoing states. The channel index m runs over all 
the transverse modes and different spin orientations. 

The operators a and b are related via an unitary trans- 
formation (n = 1, Ws) [isj : 

Ws Wo 

bsn{E)=Y^ r^,n{E)asm{E)+Y^ t' ,„n{E)aD„i{E) , 

m— 1 m— 1 

(4) 

where Ws and Wd represent the number of quantum 
channels in the source and drain leads, respectively, while 
the blocks r (size Ws x Ws) and t' (size Ws x Wd), 
describe electron reflection at the source (r) and trans- 
mission from drain to source (t') and are included in the 
scattering matrix s as (isj : 

The dimensions of s are {Ws + Wd) x {Ws + Wd)- Blocks 
t and r' in Eq. (O are related to source-to-drain trans- 
mission and reflection back to the drain, respectively. In 
the following, time dependence will be neglected, since 
we are interested to the zero frequency case. 

If we denote with | a) a many-particle (antisymmet- 
rical) state, the occupation number in the reservoir 
a in the channel m can be expressed as (Jam{E) = 
{a\^{E)aam{E)) a ■ Pauli exclusion principle does not 
allow two electrons to occupy the same spin orbital, 
therefore (JamiE) can be either or 1. In addi- 
tion, since fluctuations of the potential profile along 
the channel due to Coulomb interaction between ran- 
domly injected carriers affect the transmission of elec- 
trons, the scattering matrix elements have to depend on 
the occupation numbers of all states in both reservoirs: 
s{E) = s[asi{E),as2{E),--- ,(TDi{E),aD2{E),- ■ ■]■ Let 
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us stress the fact that, as pointed out in Ref. [13|, when- 
ever a finite channel is connected to semi-infinite leads, 
the channel can be considered as a small perturbation to 
the equilibrium regime of the contacts, and independent 
random statistics can be used for both reservoirs. 

According to Ref. current fluctuations can be 

evaluated by introducing an ensemble of many electrons 
states {|<Ti), |<T2), 1 173), • • • , I ctv)} and by weighting each 
state properly, i.e. by flnding its statistical average, de- 
noted by ( )s. Each reservoir a{a = D) is assumed to 
be at thermal equilibrium, so that its average occupancy 
can be described by the Fermi-Dirac statistics /„. As a 
consequence, the statistical average of cFam{,E) reads p^ : 

{^o.m{E))s - {{ai^{E)a^ra{E))„)s - UE) . (6) 

Neglecting correlations between the occupation numbers 
of the same quantum channel at different energies, or 
between different channels at the same energy, we ob- 
tain 

(a,™ {E)ap,, {E'))s = U {E)fp [E') (7) 

for a ^ P ov m ^ n ov E ^ E' . Including Eq. ([6]) in 
Eq. d?]) and exploiting the relation (Tam{EY = (Tam{E), 
the average of the product of two occupation numbers 
can be expressed as: 

{E)apn{E')),^ f^{E)ffi{E') + 5{E-E')5o.p5mn 
[UE)~UE)U{E')] , (8) 

where 5{E ~ E'), dap, Smn are Kronecker delta functions. 

In order to compute the average current along the 
channel and the power spectral density of the current 
fluctuations, we need to write the expectation values of 
the products of two and four operators [l^ : 

{ai^{E)apn{E'))„ = S{E - E')SapS„,nCTcm{E) ; (9) 

{ai^iE)apn{E')al,{E")asi{E"')),^ 

5{E-E"')5{E'~E")5as5ral5p^5nk(To.m{E) [1 - CT^k{E")] 

+ 6{E-E')5{E"-E"')5ap5n,n5^s5kicJam{E)a^k{E")m 

where the first contribution in Eq. (|10l) refers to exchange 
pairing (a — 6, P = ^, m = I, n = k), while the second 
to normal pairing {a — (3,^ — S,m = n, k = l) [isj . 
For the sake of simplicity, in the following we denote the 
expectation (( )„)s as ( ). 

By means of Eqs. and © the average current reads: 

•' \neS keD ) 

{neS keD ) 

(11) 



where [tj^.^^ = [fftj^^ if a = and [t'h']^^ if a = 
D (l,p € a). The unitarity of the matrix s has also 
been exploited, from which the relation r^r -I- t^t = 1 
follows. It is easy to show that for a non- interacting 
channel, i.e. when occupancy of injected states does not 
affect transmission and reflection probabilities, Eq. (fTTj) 
reduces to the two-terminal Landauer's formula 19]. 

In general, we can observe that for an interacting chan- 
nel Eq. pT|) provides a different result with respect to 
Landauer's formula, because fluctuation of transmission 
probabilities induced by random injection in the device, 
is responsible for rectiflcation of the current. The ef- 
fect is often very small, but not always [13]. However, 
it cannot be captured by Landauer's formula, as other 
many-particle processes affecting device transport prop- 
erties [illll. 

The mean squared current reads: 

{^^{lyjdEjdE' y: {{NLimLiE')) 

- {NLimsniE')) - {NsjmLiE')) 

+ {NsrnimSn(E'))} 

= (12) 

This expression consists of four terms, related to states 
at the source contacts, that can be evaluated by means of 
Eqs. (lU and ([T0| : the first one (i^++) represents the cor- 
relation of fluctuations in two ingoing streams, the second 

and the third ones {F^ , F |_) describe the correlations 

of the fluctuations of the ingoing and outgoing streams, 

the fourth one {F ) refers to two outgoing streams. 

The flrst term can be expressed as: 

E++ = {lYJdE jdE'Y. {<ys.n{E)asn{E'))s , (13) 

m , 71 G S 

since {(tI^{E))s = {asm{E))s = fs{E) Vm e S. Cor- 
relations between ingoing states are established through 
the statistical expectation values of each couple of occu- 
pancies of states injected from the source. 
The second contribution F-\ reads: 

\yndeS 

asi{E')),+Y. Y.^[t{E')]^,,^^osm{E)onk(E'))\ , 

meS keD ' ) 

(14) 

since (J^i{E) — aai{E) ^1 £ a {a — S,D), due to the 
Pauli exclusion principle. In Eq. (fT4|) correlations be- 
tween ingoing and outgoing states are obtained by sum- 
ming on each statistical average of the product of two 
occupation numbers of injected states, weighted with the 
reflection (1 — [t(£")]g^^ — [r^r{E')\ or transmission 

probability ( [t(^^')] /j.^j.) ol outgoing channels. 
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By exploiting the anticommutation relations of the 
fermionic operators a, it is simple to demonstrate that 
the third term F ^ is identical to Indeed: 



{Nsr..imUE')) = {N+^iE')N^„.iE)) . (15) 



Finally, the fourth term F reads: 



trum can be expressed as: 



a=S.D lea 



a=S.D 



l^p 



2 (^) E([t'^r],^ [vhX,-DU.S,). 

keDpes 

\ies keD / 

leSkeD 

-2 [dE IdE'Y. {[t{E)\^,^^^asi{E)asp{E')), 

•' •' i,pes 

+ (D7^/ E ME)asp{E')), . (16) 
•' •' i,pes 

Equation (1161) contains all correlations between outgoing 
electron states in the source lead, where outgoing carriers 
at the source can be either reflected carriers incident from 
S or transmitted carriers injected from D. By means of 
the Eqs. ([T^ . and p^ . we find the mean squared 
current: 



a=S,D lea 

(D'^^/^E E ([^Upi^Upi-^^-o 



i^p 



-2 {^^AEfdEj: E([t''r].p nk-n>'-sp) 
•' keD pes 

■f'^ill ra5;n^S'-E [^]D;kk''Dk 



keD 



(17) 



^ ^"^ a=S.D lea 



i^p 



Finally, from Eqs. ([T]), dTTI) and (IT71) the noise power spec- 



-2 (Q [dEj: E[t'^r],^ [rtt']^,az,.asp). 
^ ^"^ keDpes 

ERs;™'^S„-EraD;fefe^Cfc • 

I Vnes fcez? /J 

(18) 

Equation ([TS]) is the main theoretical result of this work: 
the power spectral density of the noise current is ex- 
pressed in terms of transmission (t, t'), reflection (r) am- 
plitude matrices, and properties of the leads, such as ran- 
dom occupation numbers of injected states. Let us point 
out that, although our derivation starts from Eq. ([2|), 
which is valid only far from the mesoscopic interacting 
sample, Eq. ([T5)) allows to take into account both Pauli 
and Coulomb interactions through the dependence of t, 
t' and r on actually injected states. Let us note that 
we go beyond the Hartree approximation by considering 
different random configuration of injected electron states 
for different many-particle systems. 

There is a crucial difference with respect to Landauer- 
Biittiker's formula, since Eq. enables to consider 

fluctuations in time of the potential proflle along the 
channel induced by the electrostatic repulsion between 
randomly injected electrons from the leads. Essentially, 
for each random conflguration of injected states from 
both reservoirs, we consider a snapshot of device opera- 
tion at a different time instant. All statistical properties 
— in the limit of zero frequency — can be obtained by 
considering a sufficient ensemble of snapshots. 

Let us discuss some physical limits of interest. First, 
we consider the case of zero temperature. In such condi- 
tion the Fermi factor for populating electron states in the 
reservoirs is either or 1, and all snapshots are identical, 
so the fourth term in Eq. (|18l) disappears. In addition, 
we can remove the statistical averaging in Eq. (jlSI) and 
the first three terms lead to the following expression of 
the noise power spectrum: 

5(0) ^—7- dE{TT [th] -Tr [t^tttt]) , (19) 

where Eps and Epu are the Fermi energies of the source 
and drain contacts, respectively. Such terms can be iden- 
tified with partition noise (PN) contribution. More in de- 
tail, the first term of Eq. (ITSl) is associated to the quan- 
tum uncertainty of whether an electron injected in the 
mode / from the reservoir a is transmitted through or 
reflected by the barrier. 
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The second term of Eq. ([T8)) contains instead {I p): 



[t^t]Jttt] 



pi 



E 

k,qGD 



(20) 



Each term of the sum can be interpreted as the coupHng 
between a transmission event from channel p S into 
channel k € D and from channel I G S into channel q G 
D: such a coupling is due to time-reversed transmissions 
from k into I and from q into p. 

In the same way, the third term of Eq. (|18p contains 



(21) 



that represents the coupling between carriers transmit- 
ted from n E S into k E D and reflected from p E S 
into I E S. The second and third terms provide insights 
on exchange effects. Indeed, in such terms, contributions 
with k ^ q and I ^ n, respectively, are complex and they 
represent exchange interference effects (fourth-order in- 
terference effects) in the many-particle wave- function due 
to the quantum-mechanical impossibility to distinguish 
identical carriers [13] ■ In the Results section, we will be 
concerned with identical reservoirs, i.e. identical injected 
modes from the contacts. In this case the diagonal terms 
of the partition noise (first term and part of the third 
term in Eq. (|18p ) will be referred as on-diagonal Parti- 
tion Noise (PN ON), while the off-diagonal ones (second 
term and part of the third term in Eq. HH])) will be de- 
noted as off-diagonal contribution to the partition noise 
(PN OFF). 

Now let us assume that the number of quantum chan- 
nel in the source is smaller than the one in the drain 
(Wg < Wd) and let us consider the case of potential 
barrier wide with respect to the wavelength, so that one 
may neglect tunneling. In such a situation, the reflection 
amplitude matrix r is equal to zero for energies larger 
than the barrier maximum Ec, whereas the transmis- 
sion amplitude matrix is zero for energies smaller than 
Ec- By means of the unitarity of the scattering matrix 
s, follows t^t = Is for E > Ec, where I5 is the iden- 
tity matrix of order Ws- Due to reversal time symmetry, 
there are Ws completely opened quantum channels in the 
drain contact and Wd — Ws completely closed. In this 
situation only the fourth term in Eq. (1181) survives and 
the noise power spectral density becomes: 



^(0) 



r+oo 

/ dE[fsil-fs) + fD{l-fL 
JEc 



2e^kTWs 



[fs{Ec) + Id{Ec 



(22) 



When Eps = Epo such term obviously reduces 
to the thermal noise spectrum AkTG, where G — 
(e^ Wsfs{Ec))/{'^^) is the channel conductance at equi- 
librium. The fourth term in Eq. (jl8|) can be therefore 
identified with the Injection Noise (IN) contribution. 



Equation ([18)) describes correlations between transmit- 
ted states coming from the same reservoirs [second term 
in Eq. ([TSt ] and between transmitted and reflected states 
in the source lead (third term) , with a contribution of op- 
posite sign with respect to the first term. The negative 
sign derives from Eq. (jlOp . in which exchange pairings 
include a minus sign due to the fermionic nature of elec- 
trons. Note that Eq. (IT51) can be expressed in a symmet- 
ric form with respect to an exchange between the source 
and the drain contacts. Indeed, by exploiting the unitar- 
ity of the scattering matrix, the third term becomes: 

-(^)/'^EB[t'H.p [rh']^^an,asp)s 
^ ken pes 

- /^E Et'-'^tl.p -^^-^p)^ '(23) 
\^"-^-' ken pes 

which establishes correlations between transmitted and 
reflected states in the source and drain leads. 

Now let us consider the limit when transmission and 
reflection matrices do not depend on random occupation 
numbers of injected states, i.e. a non fluctuating poten- 
tial profile is imposed along the channel. By exploiting 
the reversal time symmetry {s — s*, so that t' = t*), 
the unitarity of the scattering matrix, Eq. (jl8p reduces 
to Landauer-Biittiker's noise formula [13]: 

^(0)=^ 1 f'^ E [t^t] -Tr [tttttt] + T„) 

- [ dE 5] r„ 

a=S,D 

-2jdE (Tr [ttt] - Tr [tttttt] ) 

dE J2 (Tr[tttttt]-T„) [/„(l-/„)] I , (24) 

a=S,D J 

where T„ = Ez/pea [^^.ip Ha;^; ^'^'^ '^^^^ "^^^ 

run on the spin. Equation (j24p then reduces to: 

3(0) ^^JdE{[fsil-fs) + .fDil-fD)]Tv [thth] + 

[/5(l-/D) + /z5(l-/5)](Tr [ttt] -Tr [t^tttt])} . 

(25) 



Let us note that Eq. (fT9| can be recovered as well from 
Eq. (|25|) . Indeed at zero temperature the stochastic in- 
jection vanishes since random statistics coincides to the 
Fermi factor. In the same way, Eq. (|22p might be derived 
from Eq. (|25l). since in this case noise is only due to the 
thermionic emission contribution and fluctuations of the 
potential profile do not play any role in noise. 
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III. COMPUTATIONAL METHODOLOGY AND 
QUANTITATIVE ANALYSIS 

In order to properly include the effect of Coulomb 
interaction, we self-consistently solve the 3D Poisson 
equation, coupled with the Schrodinger equation with 
open boundary conditions, within the NEGF formalism, 
which has been implemented in our in-housc open source 
simulator NanoTCAD ViDES For what concerns 

the boundary conditions of Poisson equations, Dirichlet 
boundary conditions are imposed in correspondence of 
the metal gates, whereas null Neumann boundary condi- 
tions are applied on the ungated surfaces of the 3D sim- 
ulation domain. In particular the 3D Poisson equation 
reads 



eV(/.(f) =-(p(rO+p/„(r)) 



(26) 



where cf) is the electrostatic potential, pjix is the fixed 
charge which accounts for ionized impurities in the doped 
regions, and p is the charge density per unit volume 



r+oo 

[E) 

J Ei o r-i „ ^ ^, 



a=S,D nea 



e dEj2 ^i505„„(f,S)[l-CT„„(S)],(27) 



' a=S,D nea 

where Ei is the mid-gap potential, DOSan{r, E) is the 
local density of states associated to channel n injected 
from contact a and r is the 3D spatial coordinate. 

From a computational point of view, modeling of the 
stochastic injection of electrons from the reservoirs has 
been performed by means of statistical simulations tak- 
ing into account an ensemble of many electron states, i.e. 
an ensemble of random configurations of injected elec- 
tron states, from both contacts. In particular, the whole 
energy range of integration (Eqs. and (P?])) has been 
uniformly discretized with energy step AE. Then, in or- 
der to obtain a random injection configuration, a random 
number r uniformly distributed between and 1 has been 
extracted for each electron state represented by energy E, 
reservoir a and quantum channel n pij . More in detail, 
the state is occupied if r is smaller than the Fermi-Dirac 
factor, i.e. asn{E) [aDniE]] is 1 if r < fsiE)[fDiE)], 
and otherwise. 

The random injection configuration generated in this 
way has been then inserted in Eq. (|27l) and self-consistent 
solution of Eqs. ([26|) and ([27| and the Schrodinger equa- 
tions has been performed. Once convergence has been 
reached, the transmission (t,t') and refiection (r) ma- 
trices are computed. The procedure is repeated several 
times in order to gather data from a reasonable ensem- 
ble. In our case, we have verified that an ensemble of 
500 random configurations represents a good trade-off 
between computational cost and accuracy. Finally, the 
power spectral density S{0) has been extracted by means 
of Eq. ^ 
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FIG. 2: (Color online) a) Noise power spectral density 5(0) 
obtained from Eq. (|18p for a given potential as a function of 
current sample number for four different energy steps, b) Rel- 
ative deviation of S{0) with respect to Landauer-Biittiker's 
limit ([25} . The simulated structure is the SNW-FET shown 
in Fig. El 



In the following, we will refer to self-consistent Monte 
Carlo simulations (SC-MC), when statistical simulations 
using the procedure described above, i.e. inserting ran- 
dom occupations (Jsn{E) and <TDn{E) in Eq. (I27|) . are 
performed. Instead we will refer to self-consistent (SC) 
simulations when the Poisson-Schrodinger equations are 
solved considering fs and fo in Eq. (P7| . SC-MC sim- 
ulations of randomly injected electrons allow to consider 
both the effect of Pauli and Coulomb interaction on noise. 

From a numerical point of view, particular attention 
has to be posed on the choice of the energy step AE. 
In Fig. [5] the noise power spectrum computed by keep- 
ing fixed the potential profile along the channel and per- 
forming statistical Monte Carlo simulations of randomly 
injected electrons is shown for four energy steps. As al- 
ready proved in Eq. (IMl) . the convergence to Landauer- 
Biittiker's limit is ensured for all the considered energy 
steps: as can be seen, AE = 5 x lO"** eV provides faster 
convergence as compared to the other values with a rel- 
ative error close to 0.16%. 

Let us point out that the NEGF formalism computes 
directly the total Green's function G of the channel and 
the broadening function of the source {Ts) and drain 
(Fc) leads, rather than the scattering matrix s, that 
relates the outgoing waves amplitudes to the incoming 
waves amplitudes at different reservoirs. In order to ob- 
tain the matrix s, we have exploited the Fisher-Lee re- 
lation [2^ , which expresses the elements of the s-matrix 
in terms of the Green's function G and transverse mode 
eigenfunctions (see Appendix J. 
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FIG. 3: (Color online) 3-D structures and transversal cross 
sections of the simulated CNT (top) and SNW-FETs (bot- 
tom). 



IV. RESULTS 

The approach described in the previous section has 
been used to study the behavior of shot noise in quasi- 
ID channel of CNT-FETs and SNW-FETs with identical 
reservoirs (Fig. (H]). We consider a (13,0) CNT embedded 
in Si02 with oxide thickness equal to 1 nm, an undoped 
channel of 10 nm and n-doped CNT extensions 10 nm 
long, with a molar fraction / = 5 x 10"'^. The SNW- 
FET has an oxide thickness {tox) equal to 1 nm and the 
channel length (L) is 10 nm. The channel is undoped 
and the source and drain extensions (10 nm long) are 
doped with Nd = 10^*^ cm~^. The device cross section 
is 4x4 nm^. 

From a numerical point of view, a p^-orbital tight- 
binding Hamiltonian has been assumed for CNTs [26, 27], 
whereas an effective mass approximation has been con- 
sidered for SNWs [H, [2I] by means of an adiabatic de- 
coupling in a set of two-dimensional equations in the 
transversal plane and in a set of one-dimensional equa- 
tions in the longitudinal direction for each ID subband. 
For both devices, we have developed a quantum ballis- 
tic transport model with semi-infinite extensions at their 
ends. A mode space approach has been adopted, since 
only the lowest subbands take part to transport. In par- 
ticular, we have verified that four modes are enough to 
compute the mean current both in the ohmic and satu- 
ration regions. All calculations have been performed at 
room temperature (T = 300 K). 

Let us focus our attention on the Fano factor F, de- 
fined as the ratio of the actual noise power spectrum 
S{0) to the full shot noise 2q{I). In Figs. H and [S] the 
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FIG. 4: (Color online) Contributions to the Fano factor in a 
CNT-FET of the on-diagonal and off-diagonal partition noise 
and of the injection noise (respectively on-diagonal and off- 
diagonal part of the first three terms, and fourth term in 
Eq. (|18|l ) as a function of the gate overdrive Vas — Vth for a 
drain-to-source bias Vds = 0.5 V. a) The on-diagonal parti- 
tion (PN ON, solid circles), the injection (IN, open triangles 
up) and the full noise (open circles) computed by means of 
SC-MC simulations are shown. The Fano factor computed 
by exploiting Landauer-Biittiker's formula (|25[1 and SC sim- 
ulations (solid triangles down) is also shown, b) Off-diagonal 
partition noise contribution (PN OFF) to F due to correla- 
tion between transmitted states and between transmitted and 
reflected states. 



contributions to F of partition noise (first three terms in 
Eq. (|18|)) and injection noise (fourth term in Eq. (fT8|)) are 
shown, as a function of the gate overdrive Vgs ~ Vth for 
a drain-to-source bias Vds = 0.5 V for CNT-FETs and 
SNW-FETs, respectively: results have been obtained by 
means of SC-MC simulations. The threshold voltage Vth 
at Vds = 0.5 V is 0.43 V for the CNT-FET and 0.13 V for 
the SNW-FET. In particular, Figs. Ha andEji refer to the 
on-diagonal contribution to the partition noise (solid cir- 
cles), to the injection noise (open triangles up) and to the 
complete Fano factor (open circles) obtained by means 
of Eq. p8|) . i.e. Pauli and Coulomb interactions simul- 
taneously considered. We present also the Fano factor 
(solid triangles down) computed by applying Eq. (f25|) on 
the self-consistent potential profile, i.e. when only Pauli 
exclusion principle is included. In Figs. |4}d and [5)d we 
show the contribution of the off-diagonal partition noise 
to F, which provides a measure of mode-mixing and of 
exchange interference effects. 

As can be seen in Fig. |4^ and Fig. [5^, in the sub- 
threshold regime {Vgs ~ Vth < -0.2 V, (/) < 10^^ A) 
the Poissonian noise for a non-degenerate injection is 
recovered, since electron-electron interactions are neg- 
ligible due to the very small amount of mobile charge 
in the channel. In the strong inversion regime instead 
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FIG. 5: (Color online) Contributions to the Fano factor in a 
SNW-FET of the on-diagonal and off-diagonal partition noise 
and of the injection noise, obtained for Vds = 0.5 V, as a 
function of the gate overdrive Vas — Vth in a SNW-FET. In 
a) the on-diagonal partition, the injection and the full noise 
computed by means of SC-MC simulations (both Pauli and 
Coulomb interactions taken into account) are shown together 
with results obtained by means of Eq. (I25p . b) Off-diagonal 
partition noise due to correlation between transmitted states 
and between transmitted and reflected states. 



{Vgs -Vth>0 V, (/) > 10"'' A), noise is greatly sup- 
pressed with respect to the full shot value. In particular 
for a SNW-FET, at Vgs - Vth ~ 0.4 V ((/) « 2.4 x 
10~^ A), combined Pauli and Coulomb interactions sup- 
press shot noise down to 22% of the full shot noise value, 
while for CNT-FET the Fano factor is equal to 0.27 at 
Vgs - Vth w 0.3 V ((/) « 1.4 x 10"^ A). This is due 
to the fact that as soon as an electron is injected, the 
barrier height along the channel increases, leading to a 
reduced transmission probability for other electrons. 

As shown in Fig. |4^, the dominant noise source in bal- 
listic CNT-FETs is the on-diagonal partition noise and 
the noise due to the intrinsic thermal agitations of charge 
carriers in the contacts (injection noise), which is at most 
the 36 % of the partition noise {Vgs — Vth ~ —0.1 V). 
Nearly identical results are shown for SNW-FETs, with 
the exception of a stronger contribution given by the in- 
jection noise, up to the 86 % of the on-diagonal partition 
term {Vgs — Vth ~ —0.2 V). Moreover, the behavior of 
the two noise components, as a function of Vgs — Vth, is 
very similar for both CNT- and SNW-FETs: F tends to 
1 in the subthreshold regime, while in strong inversion 
regime shot noise is strongly suppressed. 

Let us stress that an SC-MC simulation exploiting 
Eq. (|18p is mandatory for a quantitative evaluation of 
noise. Indeed, by only considering Pauli exclusion prin- 
ciple through formula (1251) , one would have overestimated 
shot noise by 180 % for SNW-FET {Vgs - Vh ~ 0.4 V) 
and by 70 % for CNT-FET {Vgs ~ Vh ~ 0.3 V) [13, [11]. 

It is interesting to observe that the off-diagonal con- 



tribution to partition noise, due to exchange correlations 
between transmitted states and between transmitted and 
reflected states, has a strong dependence on the height 
of the potential profile along the channel (variation of 
5 orders of magnitude for CNT-FETs) and is negligible 
for quasi one-dimensional FETs. In particular, for CNT- 
FETs such term is at most 5 orders of magnitude smaller 
than the on-diagonal partition noise or injection noise in 
the strong inversion regime {Vgs — Vth ~ 0.3 V), while 
in the subthreshold regime its magnitude still reduces 
(about 10-" for Vgs - Vh ~ -0.4 V). For SNW-FETs 
we have obtained similar results: the off-diagonal parti- 
tion noise is indeed at most 5 orders of magnitude smaller 
than the other two contributions. 

In such conditions, transmission occurs only along sep- 
arate quantum channels and an uncoupled mode ap- 
proach is also accurate. Indeed, off-diagonal partition 
noise provides an interesting information on the strength 
of the mode-coupling which, as already seen, is very 
small. In particular, neglecting this term, results ob- 
tained from Eq. can be recovered as well. 

In the previous discussion, carriers from different quan- 
tum channels do not interfere. However, since we deal 
with a many indistinguishable particle system, such ef- 
fects can come into play. To this purpose, we investigate 
in more detail two examples in which exchange pairings, 
that include also exchange interference effects, give a non 
negligible contribution to drain current noise. In the past 
exchange interference effects have been already predicted 
for example in ballistic conductor with an elastic scatter- 
ing center in the channel [s^, in diffusive four-terminal 
conductors of arbitrary shape 3l| and in quantum dot 



in the quantum Hall regime |32| . connected to two leads 
via quantum point contacts. 

In the first case we discuss, mode-mixing does not ap- 
pear, i.e. the non-diagonal elements of the matrices t^t 
and t'^r are negligible with respect to the diagonal ones. 
Since the off-diagonal partition noise is negligible and 
since in the third term in Eq. (|18p only contributions 
with indices I = n — k — p survive, exchange interference 
effects do not contribute to electrical noise. We consider 
a CNT-FET at low bias condition: Vds = 50 mV. In 
Fig. [6^ the on-diagonal partition noise, the injection noise 
and correlations due to the off-diagonal partition noise, 
evaluated performing statistical SC-MC simulations, are 
shown. In this case, on-diagonal correlations between 
transmitted and reflected states in the source lead (in 
the same quantum channel) extremely affect noise. In- 
deed, at the energies at which reflection events in the 
source lead are allowed, also electrons coming from D 
can be transmitted into the injecting contact S, since 
the corresponding energy states in D are occupied and 
the barrier height is small. Instead the exchange correla- 
tions represented by the off-diagonal partition noise are 
negligible, since they are at least 5 order of magnitude 
smaller than the other three terms in Eq. ((T5)) . Note 
that the noise enhancement obtained both in the inver- 
sion and subthreshold regimes is due to the fact that at 
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FIG. 6: (Color online) a): Contributions to the Fano fac- 
tor F by the on-diagonal partition noise (solid circles), and 
the injection noise (open triangles up) as a function of the 
gate overdrive Vas — Vth, for a drain-to-source bias Vds = 
50 mV. The simulated device is a CNT-FET. The full noise 
computed by means of SC-MC simulations (open circles, both 
Pauli and Coulomb interactions taken into account) and ap- 
plying Eq. (|25p (solid triangles down, only Pauli exclusion 
considered) is also shown, b): Contributions to F by the on- 
diagonal and off-diagonal partition noise and by the injection 
noise (exploiting Eq. (|18l) ) as a function of the gate overdrive 
for a CNT-FET with a vacancy in a site at the center of the 
channel. The drain-to-source bias is 0.5 V. c): Self-consistent 
midgap potential obtained by using the Fermi statistics for 
a gate voltage Vas = 0.7 V and a bias Vds = 0.5 V. Z is 
the transport direction along the channel, X is a transversal 
direction. The simulated device is the same of b) . 



low bias the current (/) becomes small, while the noise 
power spectrum 5(0) tends to a finite value, because of 
the thermal noise contribution. 

Let now consider the situation in which modes are cou- 
pled and exchange interference effects, through the ofF- 
diagonal partition noise, contribute to drain current fluc- 
tuations. We consider the interesting case in which a va- 
cancy, i.e. a missing carbon atom, is placed at the center 
of the channel of a (13,0) CNT-FET. From a numerical 



point of view, this defect can be modelled by introduc- 
ing a strong repulsive potential (i.e. -1-8 eV, much larger 
than the energy gap of a (13,0) CNT: Egap ~ 0.75 eV) in 
correspondence of such site, thus acting as a barrier for 
transmission in the middle of the channel (Fig. |6};) . 

In Fig. |6}d the three noise sources in Eq. ([T8| (on- and 
off-diagonal partition noise, injection noise) are plotted 
as a function of the gate voltage Vgs in the above thresh- 
old regime for Vds = 0.5 V, along with the full Fano fac- 
tor computed performing SC and SC-MC simulations. 
Remarkably, in this case a mode space approach taking 
into account all modes (i.e. 13) is mandatory in order 
to reproduce all correlation effects on noise. As can be 
seen, off-diagonal exchange correlations gives rise to a not 
negligible correction to the Fano factor (w 4 % of the full 
Fano factor at Vgs = 0.8 V). We observe that such cor- 
relations are only established between transmitted elec- 
trons states (second term in Eq. (fT8)) ). while correlations 
between reflected and transmitted electron states (third 
term in Eq. (ITS)) ) are negligible since almost all elec- 
trons injected from the receiving contact D are reflected 
back because of the high bias condition. In this paper 
we have assumed phase-coherent quantum transport at 
room temperature. Our tools cannot include electron- 
phonon interaction, that a room temperature may play 
a role even in nanoscale devices. Ref. has consid- 
ered the effect of electron-phonon scattering and has ne- 
glected Coulomb interaction: they find that electron- 
phonon scattering increase shot noise in the above thresh- 
old regime, due to the broadening of the energy range of 
electron states contributing to transport. 



V. CONCLUSION 

We have developed a novel and general approach to 
study shot noise in ballistic quasi one-dimensional CNT- 
FETs and SNW-FETs. By means of a statistical ap- 
proach within the second quantization formalism, we 
have shown that the Landauer-Biittiker noise formula 
(Eq. (HH)) can be generalized to include also Coulomb 
repulsion among electrons. This point is crucial, since we 
have verified that by only using Landauer-Biittiker noise 
formula, i.e. considering only Pauli exclusion principle, 
one can overestimate shot noise by as much as 180 %. 

From a computational point of view, we have quanti- 
tatively evaluated shot noise in CNT-FETs and SNW- 
FETs by self-consistently solving the electrostatics and 
the transport equations within the NEGF formalism, for 
a large ensemble of snapshots of device operation, each 
corresponding to a different configuration of the occupa- 
tion of injected states. 

Furthermore, with our approach we are able to ob- 
serve a rectification of the DC characteristics due to fluc- 
tuations of the channel potential, and to identify and 
evaluate quantitatively the different contributions to shot 
noise. We are also able to consider the exchange interfer- 
ence effects, which are often negligible but can be measur- 
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able when a defect, introducing significant mode mixing, 
is inserted in the channel. 
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Appendix 

Let us consider a 2D channel of length L and de- 
note with X and y the longitudinal direction and the 
transverse one, respectively. If the interface between 
the lead S (D) and the conductor is defined by a;s = 
{xd = 0), Gds iyoWs) = {xd ^ 0, yo; xs = 0, ys) 
represents the wavefunction at [xd — 0,?/_d) due to an 
excitation at {xs = 0,ys)- In real space the Fisher-Lee 
relation reads: 



ihy/VnV„ 



Xm (ys) 



dyo I dys Xnivo) Gds {vdWs) 

(A.l) 



where n is a mode outgoing at lead D with velocity i;„, 
TO is a mode incoming at lead S with velocity Vm and a 
is the lattice constant along the x direction. In the k- 
representation, for a conductor of uniform cross-section, 
we can exploit a mode representation in the transverse 
direction and a plane wave representation in the longitu- 
dinal direction and (jA.l[) becomes: 



ihy/VnV„ 



Gds in,m) 



(A.2) 



where 0/35(71, m) — G^s (n, fc„; to, and fc„ is the 
longitudinal wavevector of the transverse mode n. Let 
us assume both leads to be identical and denote with 



s — 



Ts-.ii 




2Nx2N 



D — 



hv(k^) 





r_D:22 



2Nx2N 



Wn^m e S and 



Vn, TO e 13. 



where r5;ii(n,TO) 

Td:22 {n,m) = (5„ 

Generalization to a CNT-FET structure is straight- 
forward. Let us indicate with Nc and Nm the number 
of carbon atoms rings and the number of modes prop- 
agating along the channel, respectively. Since the cou- 
pling between the identical reservoirs and the channel 
involve only the end-rings of the channel, Ts and 
are {NmNc) x {NmNc) diagonal matrix and the only 
non-zero blocks are the first one and the latter one, re- 
spectively: 

hv{kn) 

Ts;ii{n,m) ^ dnm Vn, TO = 1, • • • ,iVM 



^D-N,N, {n, to) = S„ 



' L 

hv{kn) 



Vn,TO = l,--- ,iVfl,i(A.3) 



By exploiting Eqs. (jA.2p and (jA.3[) we can find the trans- 
mission (t) and reflection (r) amplitude matrix: 

t„m = i\J'TD:NcNc {n, n)GNci {n, rn) \JTs-ii {m,m) 



rnjn^-Snm + i\ Ts;ii (n, n)Gii {n,m)JTs;ii {m,m) 



(A.4) 



Since at zero magnetic field t' — t*, relations (jA.4p is 
all we need to compute the power spectral density ([T]) 
from Eq. p^ . A similar procedure has been adopted for 
SNW-FETs where, from a computational point of view, 
the channel has been discretized in a sequence of slices 
in the longitudinal direction. In this case Eqs. in (|A.4p 
are obtained as well, but replacing the number of rings 
with the number of slices. 
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